Effective interaction for p/-shell nuclei 
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An effective interaction is derived for use in the fuil pf basis. Starting from a realistic G-matrix 
interaction, 195 two-body matrix elements and 4 single-particle energies are determined by fitting 
to 699 energy data in the mass range 47 to 66. The derived interaction successfully describes various 
structures of p/-shell nuclei. As examples, systematics of the energies of the first 2 + states in the 
Ca, Ti, Cr, Fe, and Ni isotope chains and energy levels of 56 > 57 > 58 Ni are presented. The appearance 
of a new magic number 34 is seen. 

PACS numbers: 21.60.Cs, 21.30.Fe, 27.40.+Z, 27.50,+e 



The nuclear shell model has been very successful in 
our understanding of nuclear structure: once a suitable 
effective interaction is found, the shell model can repro- 
duce/predict various observables accurately and system- 
atically. For light nuclei, there are several "standard" 
effective interactions such as the Cohen-Kurath and 
the USD H interactions for the p and sd shells, respec- 
tively. On the other hand, in the next major shell, i.e., 
in the p/-shell, a unified effective interaction for all nu- 
clei in this region has not been available. The p/-shell is 
quite important for a variety of problems in nuclear struc- 
ture, such as the stability/softness of the magic number 
28, and nuclear astrophysics, such as electron capture in 
supernovae explosions. Thus, a sound, systematic and 
precise description of the p/-shell nuclei is urgent and 
important. In this Letter, we present a unified effective 
interaction which one can apply to the shell-model de- 
scription of the entire p/-shell. 

The spin-orbit splitting gives rise to a sizable energy 
gap in the p/-shell between the /V/2 orbit and the other 
orbits (P3/2, Pi/2j /5/2)j producing the N or Z=28 magic 
number. However, the excitations across the gap are im- 
portant for ground and excited state properties of many 
p/-shell nuclei. It is intriguing to understand how this 
magic number persists or fades away in various situa- 
tions. We shall discuss this point in this Letter, and 
the word "cross-shell" refers to the N or Z=28 shell gap 
hereafter. Because of such cross-shell mixing, shell-model 
calculations including all p/-shell configurations are nec- 
essary, and the predictive power obtained with a unified 
effective interaction in the full p/-shell space is very im- 
portant. The full p/-shell calculation, however, leads to 
the diagonalization of huge Hamiltonian matrices with 
dimensions of up to 2 billion. The extreme difficulty of 
dealing with such large matrices is the major reason why 
a unified effective interaction has been missing for the 
pf shell. The Monte Carlo Shell Model (MCSM) intro- 



duced recently || ^ has changed the situation by mak- 
ing such calculations feasible over the entire region of the 
p/-shell. Also conventional shell-model calculations have 
advanced. 

The effective interaction can in principle be derived 
from the free nucleon-nucleon interaction. In fact such 
microscopic interactions have been proposed for the pf- 
shell H |^| with certain success particularly in the be- 
ginning of the shell. These interactions, however, fail in 
cases of many valence nucleons, e.g., 48 Ca || and 56 Ni. 

These microscopic interactions can be modified em- 
pirically so as to better reproduce experimental data. 
The monopole-modified interaction, KB3 0, and also 
the shell-gap readjusted version, KB3G§, appear to be 
quite successful in the lower pf -shell (A < 52) . But these 
modifications turn out to be insufficient for a consistent 
description of the cross-shell properties: e.g., the 2^ level 
of 56 Ni with KB3G is predicted about 2 MeV higher than 
the observed value. 

The FPD6 interaction^] is of another type: an an- 
alytic two-body potential was assumed with parameters 
determined by a fit to the experimental data of ^4=41^49 
nuclei. It successfully describes heavier p/-shell nuclei, 
such as 56 Ni and 64 Ge |, |. There are, however, some 
defects, for instance, in the single-particle aspects of 57 Ni. 

We now turn to our new interaction. An effective in- 
teraction for the p/-shell can be specified uniquely in 
terms of interaction parameters consisting of 4 single- 
particle energies e a and 195 two-body matrix elements 
Vjriab; cd), where a, b, ■ ■ ■ denote single-particle orbits, 
and JT stand for the spin-isospin quantum numbers. 
Although e a 's include kinetic energies as well, they are 
treated as a part of the effective interaction as usual. 
We adjust the values of the interaction parameters so as 
to fit experimental binding energies and energy levels. 
We outline the fitting procedure here, while details can 
be found in Jul. For a set of N experimental energy 
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1, ■ ■ • , N), we calculate corresponding 
shell-model eigenvalues A^'s. We minimize the quantity 

X 2 = £fc=i C^exp ~ ^k) 2 by varying the values of the in- 
teraction parameters. Since this minimization is a non- 
linear process with respect to the interaction parameters, 
we solve it in an iterative way with successive variations 
of those parameters followed by diagonalizations of the 
Hamiltonian until convergence. 

Experimental energies used for the fit are limited to 
those of ground and low-lying states. Therefore, cer- 
tain linear combinations (LC's) of interaction parameters 
are sensitive to those data and can be well determined, 
whereas the rest of the LC's are not. We then adopt 
the so-called LC method [0, where the well-determined 
LC's are separated from the rest: starting from an initial 
interaction, well-determined LC's are varied by the fit, 
while the other LC's are kept unchanged (fixed to the 
values given by the initial interaction). 

In order to obtain shell-model energies, both the con- 
ventional and MCSM calculations are used. Since we 
are dealing with global features of the low-lying spec- 
tra for essentially all p/-shell nuclei, we use a simpli- 
fied version of MCSM: we search a few (typically three) 
most important basis states (deformed Slater determi- 
nants) for each spin-parity, and diagonalize the Hamilto- 
nian matrix in the few-dimensional basis approximation 
(FDA). The energy eigenvalues are improved by an em- 
pirical correction formula with parameters fixed by cases 
where more accurate results are available. This is called 
hereafter few-dimensional approximation with empirical 
corrections (FDA*), which actually yields a reasonable 
estimate of the energy eigenvalues with much shorter 
computer time. 

In the selection of experimental data, in order to elim- 
inate intruder states from outside the present model 
space, we consider nuclei of A > 47 and Z < 32. As 
a result 699 data of binding and excitation energies (490 
yrast, 198 yrare and 11 higher states) were taken from 

87 nuclei : 47 " 51 Ca, 47 - 52 Sc, 47-52^ 47-53,55y ; 48-56^ 
50-58 Mn; 52-60^ 54-61 C ^ 56-66^ 58-63 CU; 60-64^ 

62,64,65Q a anc j 64,65Q e> ^y e assume an empirical mass 
dependence A~ 3 of the two-body matrix elements sim- 
ilarly to the USD interaction^. We start from the real- 
istic G-matrix interaction with core-polarization correc- 
tions based on the Bonn-C potential which is simply 
denoted G hereafter. 70 well-determined LC's are varied 
in the fitting procedure, and a new interaction, GXPF1, 
was obtained with a rms error 168 keV within FDA*. 

The first 2 + energy level of even-even nuclei are a good 
systematic measure of the structure. The left panel of 
Fig. H shows them for the Ca, Ti, Cr, Fe, and Ni isotopes. 
The lightest nucleus in each isotope chain corresponds to 
N=Z because of the mirror symmetry. The energies are 
computed by FDA* and by exact or nearly exact conven- 
tional shell- model calculations by the code MSHELL[p"2||. 
The former gives a reasonable approximation to the lat- 
ter, while the differences are up to 0.2 MeV. The overall 
description of the 2 + levels is quite successful throughout 
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FIG. 1: (Left) First 2 + levels as a function of the neutron 
number N. Experimental data are shown by filled circles |l| 
and open triangles [Q. Solid lines show results of conven- 
tional calculations: the maximum number of nucleons excited 
from f 7/2 to p 3/2 , Pl/2 or f 5/2 is 5 for 56 Fe and 58 > 6 °. 62 Ni, 6 
for 52 - 54 Fe and 56 Ni, and 7 for 58 ' 60 Fe, whereas the others are 
exact. Dashed lines imply FDA* results. (Right) Effective 
single-particle energies for neutron orbits. Symbols indicate 
that the corresponding orbit is occupied by at least one nu- 
cleon in the lowest filling configuration. 



these isotope chains. In all cases, the energy jump corre- 
sponding to N — 28 shell closure is nicely reproduced. 

A basic aspect of the effective interaction is provided 
by the effective single-particle energies (ESPE)|, ||. 
The ESPE depends on the monopole part of the Hamil- 
tonian, and reflects angular-momentum-averaged effects 
of the two-body interaction for a given many-body sys- 
tem. In the right panel of Fig. [j], ESPE's of the neutron 
orbits are shown. A new A^=34 magic structure has been 
predicted in In fact a large energy gap (^4.1 MeV) 
between P1/2 and f§/2 can be seen in the ESPE for Ca 
and Ti isotopes. 

As predicted also in [^6|, due to the large attractive 
proton- neutron matrix elements, the ESPE of the 1^/5/2 
orbit comes down as 7r/ 7 / 2 orbit is occupied, reducing this 
subshcll gap. One can indeed see this change in Fig. |]: 
the A^=34 subshell becomes weaker in Cr, and disappears 
in Fe and Ni. This means that the A^=34 magic number 
arises only in neutron-rich Ca and Ti isotopes, and is 
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FIG. 2: 



Energy levels of 56 > 57 . 58 Ni. Experimental data are taken from ref.( 
only for yrast states for 57 ' 58 Ni. 



Above 4 MeV, experimental levels are shown 



gone in stable nuclei. On the other hand, the ESPE's 
of Ni isotopes in Fig. |l| indicate that the f$/ 2 , ^3/2,1/2 
orbits are degenerate to a good extent, forming a subshell 
which resembles a degenerate pseudo sd- shell. This fact 
may be relevant in constructing some algebraic models, 
for instance, a pseudo SU(4) of IBM-4 |Tj. 

In the Ca isotopes, a prominent peak in the calculated 
2 + excitation energy can be seen at 7V=34, which is as 
high as that of the doubly magic nucleus 48 Ca. This is ex- 
actly due to the iV=34 subshell closure discussed above. 
Note that this calculation is a result of the diagonaliza- 
tion of the Hamiltonian, whereas the ESPE reflects only 
its monopole part. The subshell gap at N—32 is much 
smaller (~1.6 MeV) between p 3 / 2 and P\/ 2 , and shows 
less pronounced effect in a consistent manner with ex- 
periment. The gap at N = 34 is not large with the 
FPD6 interaction as can be inferred from the 2 + state 
systematics shown in Fig. 4 of [Il8f . Thus the experimen- 
tal energy of the 2 + state in 5i nTa will be an important 
test of the pj-shell Hamiltonians. 

In the Ni isotopes, both experimental and theoretical 
2 + excitation energies drop at 7V=34, where ESPE of 
vpi/2 and f/5/2 are almost degenerate and therefore the 
collectivity is enhanced. It is remarkable that the drastic 
change of the structure among those nuclei can be de- 
scribed by a single effective interaction. Toward the end 
of the p/-shell (iV=40), effects of g 9 / 2 seem to appear. 

The nucleus 56 Ni is a challenge for a unified descrip- 
tion of nuclei in the cross-shell region of the p/-shell. 
The stability of the (/7/2) 16 core plays important roles, 
while most of existing effective interactions fail in repro- 
ducing certain properties related to the core softness. In 
Fig. g energy levels of 56 > 57 > 58 Ni are shown. Theoretical 
levels are obtained by MCSM calculations, where typi- 
cally 13 J-compressed basis statesjQ are taken for each 
state. These MCSM calculations are more accurate than 
the FDA*. The agreement between the theory and the 
experiment is satisfactory. 



In the calculated spectrum of 56 Ni, the deformed band 
discussed in 1 19, 2?J appears as Of, 2 2 , 4^ , 6^, • • ■ . In the 
yrast band, we obtain B(E2; 0^ -> 2+)=5.5xl0 2 e 2 fm 4 , 
while in the deformed band, B(E2; 0+ -> 2+)=1.8xl0 3 
e 2 fm 4 . Here, effective charges, e p =1.23, e n =0.54, are 
taken Q. The former _B(E2) value is in agreement with 
experiment |2l]] . The probability of the (/V/2) 16 con- 
figuration in the ground state is 69%. This value is 
larger than that obtained by FPD6 (49%) B, while is 
smaller than the corresponding quantity for Ca (94% 
for GXPF1). 

In order to study the core-softness, we made a rather 
"soft" interaction (GXPF2). For 20 matrix elements of 
the type VjT(aa]bc) where a — j-jj 2 and b,c ^ a, we 
keep their G interaction values. We then carried out 
a fit for 60 best determined LC's, and came up with a 
rms error of 188 keV for 623 data within FDA*. This 
means that the fit is only slightly worse (by ~20 keV) 
than GXPF1. The probability of (77/2) 16 configuration 
in the ground state of 56 Ni is 49% suggesting its soft- 
ness. Thus the energy data included in the fit can be 
reproduced about equally well within certain allowance 
of the softness. The core property, however, is important 
for certain observables. For example, the total Gamow- 
Teller strength from the ground state is calculated as 
B{GT+)=ll.i (9.5) for GXPF1 (GXPF2). Since this 
is 13.7 for the closed shell, one sees the degree of the 
quenching which can be crucial in the electron capture 
in the supernovae explosion. From the KB3, one obtains 
10.1 @. 

In the lowest three states 3/2~, 5/2~, and l/2~ of 
57 Ni, the (/7/2) 16 core is broken similarly to 56 Ni ground 
state, and these are "single-particle" states built on top 
of the correlated 56 Ni ground state. Their relative posi- 
tions are determined mainly by the ESPE's. In pj-shell 
nuclei, in general, the cross-shell excitations occur rather 
commonly in low- lying states p3l. In 57 Ni, states above 
the three "single-particle" states contain further cross- 
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FIG. 3: Energy levels of odd-A Ni isotopes. Experimental 
data are taken from ref.J^]. Theoretical results are calcu- 
lated by using GXPF1/2 interactions, which are indicated as 
1 and 2, respectively. Shell model calculations were carried 
out by the code MSHELL(i|]. The results for 57 - 63 Ni are ob- 
tained with truncating the number of nucleons excited from 
fr/2 orbit up to six, whereas the result for 65 Ni is exact. 



shell excitations, which provide a good testing ground for 
cross-shell two-body matrix elements. A similar situation 
has been seen in 56 Ni particularly for its 4p-4h deformed 
band, while such effects are less evident in 58 Ni, with an 
exception of the state which contains a sizable amount 
of proton 2p-2h excitations 0]. In these examples np- 
nh refers to the particle-hole excitations across the gap 
on top of the correlated ground state A compre- 

hensive picture for such a wide variety of states requires 
the full-space calculations with an appropriate effective 
interaction. 

The description of odd- A nuclei is also an important 
test of the interaction. We again consider Ni isotopes as 
examples, for which energy levels of yrast l/2~, 3/2 ~, 
and 5/2~ states are shown in Figj3[ Experimentally, as 
the number of neutrons increases, both 1/2 - and 5/2 _ 
states come down relative to 3/2 ~ level, forming nearly 
degenerate states around 63 Ni. This characteristic fea- 
ture is nicely reproduced by both GXPF1/2 interactions. 
The downward slope of 5/2~ level, however, seems to be 
somewhat too steep, especially in GXPF1, which leads to 
some deviations within about 0.2 MeV, while the basic 
degeneracy around 63 Ni is clearly maintained. Since the 
results of the FDA* used in the fit are much closer to 
experimental spectra, a part of this deviation is due to 
the uncertainty in the FDA*. 

These yrast levels contain similar amount of cross-shell 
excitations within each isotope. In fact, the occupation 
number of / 7 / 2 orbit varies from 15.0 (14.6) for 57 Ni to 
15.6 (15.5) for 65 Ni by GXPF1 (GXPF2), while the dif- 
ference among the lowest 3 states within each isotope is 
less than 0.3. Therefore the core excitation appears to 
play a minor role for describing relative energies of these 
lowest 3 levels. However, it should be emphasized that 
the GXPF1/2 interactions can describe these lowest lev- 



TABLE I: Single-particle energies relative to p 3 / 2 of various 
interactions. 
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els as well as non-yrast levels for which the cross-shell 
excitation is essential. 

In usual shell model interactions, we often take the 
single-particle energies e a from experimental energy spec- 
tra of the one-particle/hole system on top of the assumed 
inert core. In the case of p/-shell, many existing effective 
interactions such as KB3 borrow them from 41 Ca. How- 
ever, our main purpose is not to describe light pj-shell 
nuclei from the beginning of the shell, but to treat cross- 
shell excitations over N or Z = 28 shell gap. Therefore, 
in the present approach, we have assumed single-particle 
energies as parameters and determined them by the fit. 

In Table | the single particle energies relative to p 3 / 2 
orbit are shown for various interactions. It is remark- 
able that the single-particle energy spacing between f? / 2 
and p 3 /2 orbits in GXPF1/2 interactions is enhanced by 
about 0.9 MeV in comparison to the energy spectrum of 
41 Ca (or KB3). This difference in single-particle ener- 
gies results from the exclusion of energy data of A < 47 
nuclei from the fit. The wave functions of nuclei near 
40 Ca contain relatively large amounts of intruder state 
admixtures, and the sd-shell should be included in their 
description. A fit in which the single-particle energies 
are fixed to their values in 41 Ca is possible, but the total 
rms error in the energies becomes larger. The ESPE of 
GXPF1 interaction for Ca isotopes becomes very close to 
that of KB3, which is known to be quite successful for 
light p/-shell nuclei, already around A ~ 45 (where the 
difference is less than 0.2 MeV). Thus the description of 
A < 47 nuclei by GXPF1/2 interactions is still reasonable 
and similar to that of KB3, except for the those nuclei 
at the very beginning of the shell. 

In Fig. |^(a), a comparison between GXPF1 and G is 
shown for the 195 two-body matrix elements. One finds 
a strong correlation. On average, the T=0 (T=l) matrix 
elements are modified to be more attractive (repulsive). 
The majority of most attractive matrix elements are T=0 
ones, and, in particular, the two most attractive ones be- 
long to T=0 f7f5 in both GXPF1 and G, where the nota- 
tion f7f5 refers to the set of matrix elements Vj^r(ab\ ab) 
with a = fr/2, b = f 5/2 - 

In Fig. H(b) the difference between GXPF1 and G is 
shown for several T—Q diagonal matrix elements. The fit 
makes them more attractive. Although the f7f7 matrix 
elements are dominated by the monopole shift, i.e., a 
J-independent correction, the corrections to cross-shell 
matrix elements f7p3 and f7f5, become larger for higher 
J's. Such J-dependences beyond the monopole shifts are 
of interest. The same feature is also found in GXPF2. 
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FIG. 4: (a) Correlation of V JT (ab; cd) between G and GXPF1. 
T=0 and T=l matrix elements are shown by open circles and 
crosses, respectively, (b) Difference of diagonal T—0 matrix 
elements between G and GXPF1 as a function of the angular 
momentum, J, coupled by two nucleons. 

Apart from these corrections, notable differences between 
G and GXPF1 are found in the T = 1 monopole-pairing 
cross-shell matrix elements. 

In summary, a new unified p/-shell effective interac- 
tion, GXPF1, has been obtained. The GXPF1 inter- 



action has monopole properties enforced by the energy 
data, and it properly handles cross-shell excitations, lead- 
ing to a successful description of structure of Ni isotopes. 
Corrections beyond the monopole shifts are important. 
The systematic behavior of the energies for the lowest 
2+ levels of even-even nuclei are in good agreement with 
experiment, suggesting that collective properties are well 
described. We investigated the possible range of the soft- 
ness of N=Z=28 core, coming up with another interac- 
tion, GXPF2. The applications to unexplored regimes 
of large neutron numbers or high excitation energy is 
of great interest. For example, the GXPF1 interaction 
demonstrates the appearance of a new magic number 
-/V=34 |l6[ ] in neutron-rich nuclei. Future experiments 
will test the predictions and provide guidance for further 
improvements in the Hamiltonian. 
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